ON SYSTEMS OF RATIONAL DIFFERENCE EQUATIONS AND 
PERIODIC TETRACHOTOMIES 
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Abstract. We study the fc*^ order system of two rational difference equations 

PkXn-k + IkUn-k 



SkXn-k + f-kVn-k ^ 

with nonnegative parameters and nonnegative initial conditions. We establish the 
existence of periodic tetrachotomy behavior which depends on the matrix 

Pk Ik 
5k e/c 

We provide some partial results for similar systems of three or more rational difference 
equations. 
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1. Introduction 

Recently several papers discussing rational systems in the plane have appeared in 
the literature. We refer particularly to [2], [3], and [5]. In [2] the authors mention 
a conjecture regarding periodic trichotomy behavior for some rational systems in the 
plane. Our goal here is to try to extend a general periodic trichotomy result presented 
in [20] to systems of rational difference equations of order greater than one. 
In [20] the author presents a broad collection of periodic trichotomy results for the k^^ 
order rational difference equation with non-negative parameters and non-negative initial 
conditions 

Xn = —T-^ , neN. 



The author focuses on two ideas. The first idea involves expanding upon the results 
in [19] so as to create a general periodic trichotomy result. This idea applies in the case 
where a = and A > 0. This results in a trichotomy which depends on a comparison 
between E\^i Pi and A. The second idea is a nontrivial observation. The author observes 
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that there is a general trichotomy result which follows immediately from the current 
literature with only some minor technical work required. The author then patches 
together results from [4j, [17J, [19j, |21], and |22j to create a second general periodic 
trichotomy result. This idea applies in the case where a > and also depends on a 
comparison between X]i=i A and A. In this way the results presented in |20) complement 
each other nicely. 

Let us look more closely at the first idea in |20j. The author first places conditions on 
the rational difference equation so that there is a nontrivial subspace of initial conditions 
where the solution behaves linearly. The author then shows that the rational difference 
equation inheirits trichotomy behavior from the associated linear difference equation. 
In this paper we will build off of this idea, with the hope of creating some basic examples 
of trichotomy behavior for systems of rational difference equations. 

2. Some Linear Algebra Review 

Recall that if we have a symmetric matrix with real coefficients then such a matrix 
must be Hermitian. Any such matrix A is diagonalizable and has decomposition UDU* 
where D is a diagonal matrix consisting of the eigenvalues of A, U is a unitary matrix, 
and U* represents the conjugate transpose of U. Furthermore we know that D has only 
real entries. Thus we have the following facts. 

Fact 1. Suppose we have a real symmetric m x m matrix A whose spectral radius is 1 
then there exists an N E N so that for every L > N we have {A^v, A^v) < {v,v) for all 
V G M™. Moreover {A^v, A^v) = {v,v) if and only if v is in the span of the eigenvectors 
of A with corresponding eigenvalues whose absolute value is 1. 

Fact 2. Suppose we have a real symmetric m x m matrix A whose spectral radius is 
greater than 1 and a vector v G M*" so that {v, Wi) ^ for all i G {I, . . . , m}, then the 
sequence {{A^v, A^v)}^^^ is unbounded. 

3. A Lemma 

Here we present a lemma in order to simplify the results presented in the following 
sections. 

Lemma 1. Let V be an inner product space, let A : V ^ V be a bounded linear operator 
with respect to the inner product norm, \ \ ■ ||. Assume that {fnl^i 'is bounded. Also 
assume that there exists M > so that || < Af| It'll for all n E N and for all v E V. 
Let lim„^oo \ \vn — ^^n-fc|| = and suppose ^"+™ — — )■ for some m G then 

lim^_^oo W'^n f^n— femll — 0. 

Proof. For any p G N and any n > km + kp the triangle inequality gives us 

ll^n t'n— fcmll ^ ll^'^'^n— fcp A^ Vji—km—kp\ \ ~l~ 11^'^ '^n—km—kp ^'''^n— fcm— A;p 1 1 ~l~ 
ll^n— fcm ■^^'^n—km—kpW ~l~ ll^n ^'^^n— fcp||' 

We want to show that given e > there exists G N so that ||fn — "yn-fcrnll < e for 
n> N. Since {vn}'^=i is bounded and A""^™ — A" — )■ 0, given e > there exists G N 
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SO that \\A''^~^"^Vn-km-kpe ~ AP^Vn-km-kpeW < f ^ > + ^Pe- For the sake of 

notation, from now on we let = p. So we have 

ll'^n '^n— femll ^ fep A^ Vn—km—kp\\~^'^~^\\'^n—km ^'''^n— fem— fep| | I'^n ^'''^n— 

for any n > km + kp. Notice that the triangle inequality, and the fact that ll^"^!! < 
M\ \v\ I for all 77. e N and for all e F gives us 

m—1 

||^'''^n-fcp ~ A^'^^Vn-km-kpW — ^ ] \ \AP'^^Vn-k{p+i) ~ AP'^^'^^Vn-k{p+i+l)\\ 

i=0 

m—1 

< \ \Vn-k{p+i) - AVn-k{p+i+l)\\-, 

i=0 

for all n > km + kp. Since hm„_^oo I kn — ^^n-fcl I = 0, given e > there exists an Ni eN 
so that \\vn - Avn-k\\ < for n > Ni. So we get 

m—1 
i=0 



for n > Ni + km + kp. Thus 

2 



'^n "^n—kmll „ ~l~ ||^n— fem ^'''^n— fcm— fep| | ~l~ ||^n -'^-^''^n- fep| | > 



for any n > Ni + km + kp. 

We use an argument which is similar to the prior argument. Again the triangle inequality, 
and the fact that | l^'^w] | < M\ \v\ \ for all n e N and for all w e gives us 

p-l p-l 

\\Vn - A^Vn-kpll < y^J\A'Vn-ki " A'^'^Vn-k-ki\ \ < M^\\Vn-ki - AVn-k-ki\\, 
i=0 i=0 

for any n > km + kp. Since lim^^oo ll^n — ^^n-fc|| = 0, given e > there exists an 
N2 eN so that \\vn — Avn-k\\ < for n > N2. So we get 



\Vn - A^Vn-kpW <MJ2 



p-l 

e 



for all 77, > A^2 + km + kp. 

Again, We use a similar argument. The triangle inequality, and the fact that < 
M\ \v\ \ for all n e N and for all v G V" gives us 

p-l 

ll'^n— fem ^'''^n— fem— fepi I ^ ^ ^ ||^ f^n— fcm— fei A Vn—k—km—ki\\ 

i=0 

p-l 



— ^ ^ ll'f^n— fem— fei AVn—k—km—ki\ 



i=0 
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for any n > km + kp. Since lim„_^oo \ \vn — ^Vn-k\\ = 0, given e > there exists an 
A^2 G N so that \\vn — Avn-k\ \ < jMp ^ — -^2- So we get 

p-i 



for all n > A^2 + km + A;p. 

So we have that given e > there exists N — Ni + N2 + km + kp^ + 1 so that 

II 11 e e e 

I \Vn - Vn-km\ l'^2"^4"^4^^' 

for any n> N. Thus lim„_>oo lkn - Vn-km\ \ =0. □ 



4. A Representation Using Vector Spaces 
Let A; e N, A; 7^ 1 and consider the k*^ order system of two rational difference equations 

_ PkXn-k + IkVn-k ^ 

with nonnegative parameters and nonnegative initial conditions. We find that it is useful 
to rewrite our system using matrix notation. We let 

Xn\ A = ( 
Vn J ' \ h ^k 



and 



where 







-Dn - n 1 



Our system then becomes 



Vn = BnAVn~k, U E N. 



5. Some General Results 

Notice that in the last section we were able to represent our A;*'^ order system of two 
rational difference equations as a recursive system with delays on a two-dimensional 
vector space. In this section we prove some general results regarding recursive systems 
with delays on vector spaces of arbitrary finite dimension. 
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Theorem 1. Consider the k*^ order recursive system on [0,oo)"^ 

where A = (aij) is a real symmetric m x m matrix with nonnegative entries aij > and 
with initial conditions in [0,00)"^. Further assume that B.^ is a real rn x m diagonal 
matrix which may depend on n and on prior terms of our solution {vn}, with all entries 
bn,ij £ [Oj 1] for all n E N. Then we have the following 

I Whenever the spectral radius of A is less than 1, then every solution converges to 
the vector. 

II Suppose the spectral radius of A is equal to 1 then we have lim„_^oo | I'^^n— ^'^n-fcl | = 0, 
where \ \ ■ \ \ represents the inner product norm. 

Proof. Let us first prove case (I). Consider the system 

Un = Aun-k, n eN. 

Suppose Vn = Un ioY n < 1. In other words suppose that the two systems have the same 
initial conditions. Then the ith entry of the vector Vn is less than or equal to the ith 
entry of the vector m„ for all n G N and for alH G {1, . . . , m}, in other words Vn,i < Un,i- 
We prove this by strong induction on n. The initial conditions provide the base case. 
Suppose the result holds for n < N. 

m m m 

VN,i — bN,ii dijVN-kJ < 0,ijVN-k,j ^ dijUN-kJ — Un,^ 
j=l j=l j=l 

since bN,ii G [0, 1] and Uij > for alH, j G {!,..., m}. Thus we have shown Vn,i < Un,i 
for all n G N. 

It is clear that u^n+b = A"-Ub. Now if the spectral radius of A is less than one it is a 
well known result that lim^^oo A^ — 0. Of course by here we mean the zero matrix. 

(°\ 

Thus in this case hm„_^oo i^n = : ■ Since Vn G [0, 00)"* for all n G N we have 

V / 

limn^ooVn =1 : • This finishes case (I). 
\ / 

Now let us consider case (II). Since A is real m x m matrix with nonnegative entries 
Qij > 0, and A also has the property that every eigenvalue A with |A| = 1 has algebraic 
multiplicity equal to geometric multiplicity, Fact 1 applies. 

Now using a similar argument as was used in case (I) we will show something slightly 
more general. Let Un,q,L = A'^^^Vn-kL, we will show that Un,qfl,i < Un,q,L,i whenever 
L > for all i G {1, . . . , m}, g G N, and for all n G N with n > kL. Now we will use 
induction on q with the case q = providing the base case. In the case g = we must 
show that Unfl,o,i < Unfi,L,i for all i G {l,...,m} and whenever L > and n > kL. 
Consider the system 

Un = Aun-k. n G N. 
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Fix an L > and fix an Ni > kL and suppose v^-^^ki = uj^f^^^L- Then the ith entry of 
the vector v^-^^ is less than or equal to the ith entry of the vector Ujvi for alH G {1, . . . , m}, 
in other words VNi,i < UNi,i- We prove that 

VNi-kL+kp,i ^ U'Ni-kL+kp,i, 

by induction on p. The case p — provides the base case. Suppose the result holds for 
p < P. Then we have 

m m 

VNi-kL+kP,i = ^ATi-fcL+fcP.ii ^ aii'L'Ari-fcL+fc(P-l),i < CiijVNi-kL+k(P-l),j 

j=l 3=1 

m 

— ^ %'"A^i-fcL+fc(P-l)J = UNi-kL+kP,i, 

since hj^^^kL+kP,ii £ [0,1] and a^j > for all i,j e {l,...,m}. Thus we have shown 
in this case that t'jv^^i < u^q^^i for all i G {1, . . . I'm.}- Notice that fjv^^j = Ujsf^fifi^i and 
'WAri,j = Uiq-^fl^L^i. Further notice that this is true for any fixed L > and fixed Ni > kL. 
Thus [/n,o,o,i ^ Unfl^L,i for all i G {1, . . . ,m} and whenever L > and n > kL. So we 
have shown the case q — 0. 

Now we use the case g = as a base case in an induction argument on q. Suppose that 
Un,q,o,i < Un,q,L,i, then siucc all entries a^j of our matrix A are non-negative and Un,q,L is 
a vector with all non-negative entries, 

m m 

Un,q+l,0,i = Clij Un,q,0,j ^ (lijUn,q,L,j = Un,q+l,L,i- 
j=l j=l 

For shorthand we define a function h : — >■ [0,oo) so that h{v) — {v,v). Fact 1 
gives us that there exists an G N so that h{A^v) < h{v) for all v G and all L > N. 
Thus fixing an L > iV and a g G N we have h{A'^Vn) < h{A^~^'^Vn-kL) < K^'^'^n-kL)- 
Since each of the subsequences {h{A'*VnkL+a)}'^=i are monotone decreasing and bounded 
below by zero, they all converge. So lim^^oo h{A'^Vn) — hi^A'^Vn-kL) = 0. By the squeeze 
theorem we get lim„^oo h{A'^Vn) — h{A^~^'^Vn-kL) = 0. We use the triangle inequality and 
we get \h{vn)-h{Avn-k)\ < \h{A^'+^Vn-k~kLi) -h{Avn-k)\ + \h{A^'+^Vn-k-kLj -h{vn)\. 
Since, lim„_>oo h{A'^Vn) — h{A^^'^Vn-kL) = 0, we may substitute first with L = Li > N 
and q = I. This shows us that lim„^.oo \h{A^''-^^Vn-k-kLi) — h{Avn-k)\ = 0. We may also 
substitute L — Li + 1 and g = to see that lim„_^oo \h{'^^^^^'^n-k-kLi) — h{'^n)\ — 0. 
This implies lim.„^oo h(vn) - h{Avn-k) = 0. 

Using the notation we have established earlier we may rewrite Avn-k as f/„^o,i and we 
may write Vn = BnUn^^i. Recall that i?„ is a diagonal matrix which may depend on n 
and on prior terms of our solution {vn}, with all entries bn,ij G [0, 1] for all n G N. Thus 
since 5„ is a diagonal matrix we have 

m 

7?,— vno ^ ^ ' ' ' ' 
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Since all entries fe„ jj G [0,1] for all n G N and all entries of ?7„,o,i are nonnegative for all 
neN we have that (6^ - l)[/^^o,i,i ^ ^r alH G {1, . . . , m}. Thus 



n— >oo ' ' ' ' 

for alH G {1, . . . , m}. Thus 

lim (1 - 6^,,Xo,i,i = 0, 
for alH G {1, . . . , m}. Since 6„_jj G [0, 1] for all n G N we get, 

bl,, - 1 < {bn,u -lf<l- bl,„ 

for all n G N. This gives us 

lim {bn,u - l)Un,0,l,i = 0, 
n— ^oo 

for all i G {1, ... , m}. This implies that 

hm Vn - AVn-k = 0. 

So 

lim \\Vn - AVn-k\ \ = 0. 
n— >-oo 

This completes the proof of the theorem. □ 

Theorem 2. Consider the k*^ order recursive system on [0, oo)"* 

v„ = BnAvn_k, neN, 

where A = (aij) is a real symmetric m x m matrix with nonnegative entries Uij > and 
with initial conditions in [0,00)™". Further assume that Bn is a real m x m diagonal 
matrix with entries hn^u = — ^k-i for all n G N. Where the vectors qij G [0, oo)'". 

If the spectral radius of A is greater than 1 . Then for some choice of initial conditions 
the solution {vn}'^=i is such that {||^'n||}^i is an unbounded sequence. 

Proof. Before we begin to prove the first case notice that if we choose initial conditions 
so that = for 77. < 1 and n 1 — k, then it is clear by a simple induction argument 
that Vn — ior n ^ 1 mod k. Thus for solutions with these initial conditions we have 
Vn = Avn-k- So here we intend to take advantage of this linearity so we will assume 
that f „ = for n < 1 and n 1 — k, and our goal will be to choose vi-k appropriately 
in order to create an unbounded solution. 

If we choose Vi^k £ [0, oo)"* so that for all the eigenvectors of ^4, wi, . . . w^, {vi-k, Wi) 7^ 
for all i G {l,...,m}, this is certainly possible since [0,00)™" is an m-dimensional 
subspace of M™. We apply Fact 2 and we notice that Hi'fcL+ill = | l^^^^'^'^i-fcl | = 

^ I^A^+'^vi-k-, A^^^vi-k), thus {||t'fcL+i||}L=i is unbounded, so {||fn||}^i is unbounded. 

□ 
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6. A General Periodic Trichotomy Result 

In this section we use the Perron-Frobenius theorem along with our work in the last 
section to demonstrate a general periodic trichotomy result. For more details regarding 
the Perron-Frobenius theorem see [TOJ chapter 8 sections 2 and 3. 

Theorem 3. Consider the k^^ order recursive system on [0, oo)™ 

Vn = BnAVn-k, G N, 

where A = {oij) is a real symmetric mx m matrix with positive entries Oij > and with 
initial conditions in [0, oo)™. Further assume that is a real m x m diagonal matrix 
with entries hn/d = — ^.i^ for all n G N. Where the vectors Qij G [0, oo)"^. Then 

this system displays the following trichotomy behavior: 

i If the spectral radius of A is less than 1 then every solution converges to the vector. 

ii If the spectral radius of A is equal to 1 then every solution converges to a solution 
of not necessarily prime period k. Furthermore in this case there exist solutions of 
prime period k. 

Hi If the spectral radius of A is greater than 1 then for some choice of initial condi- 
tions the solution {vn}^=i has the property that {\\vn\\}^=i is an unbounded sequence. 
Moreover if we consider the sequences consisting of the entries ofvn, {wn,j}^i, then 
{vn,i}'^=i is an unbounded sequence for every i G {1, . . . , m}. 

Proof. First notice that (i) follows immediately from Theorem 1. Now consider case (iii). 
From Theorem 2 we get immediately that there is some choice of initial conditions so 
that the solution {vn}^=i has the property that {ll'i'nll}^! is an unbounded sequence. 
Recall from the proof of Theorem 2 that every unbounded solution we constructed had 
the property that f „ = for n < 1 and n 1 — k. For our purposes we will choose 
an unbounded solution which has this property, thus Vn = Avn-k for our solution. 
Since {||^'n||}^i is an unbounded sequence it follows as a consequence {f„ j^j^^^ is an 
unbounded sequence for some Zi G {1, . . . , m}. So there is a subsequence {vni^,ii} which 
diverges to oo. Since A = (aij) is a real m x m matrix with positive entries aij > 
and fn^+fc = Avn^, the subsequence {vn]^+k,i} diverges to oo for all i E {I, . . . ,m}. So 
{Vn,i}'^=i is an unbounded sequence for all i G {1, . . . , m}. This concludes the proof of 
case (iii). 

To prove case (ii) we use the Perron-Frobenius theorem. The Perron-Frobenius theorem 
tells us that if ^4 = {atj) is a real m x m matrix with positive entries aij > 0, then 
there is a positive real number r called the Perron-Frobenius eigenvalue such that r is 
an eigenvalue of A and so that any other possibly complex eigenvalue A has |A| < r. 
Moreover r is a simple root of the characteristic polynomial and there is an eigenvector 
Wr associated with r having strictly positive components. Now combining this with the 
fact that the spectral radius is 1 we get that r = 1 and every other eigenvalue A has 
|A| < 1. Also we know that r is a simple root of the characteristic polynomial so r has 
algebraic multiplicity equal to 1. So it must be true that every eigenvalue A with |A| = 1 
has algebraic multiplicity equal to geometric multiplicity. Thus Fact 1 applies in this 
case. 
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The next part of this argument wiU be almost identical to the argument presented in 
Theorem 1. Let Un,q,L = ^'^^^Vn~kL, we will show that Un,qfi,i < Un,q,L,i whenever L > 
for alH e {1, . . . , m}, g e N, and for all n e N with n > kL. Now we will use induction 
on q with the case q — providing the base case. In the case g = we must show that 
Un,o,o,i < Unfi,L,i for all i G {1, . . . , m} and whenever L > and n > kL. Consider the 
system 

Un = Aun-k, n eN. 

Fix an L > and fix an A^i > kL and suppose vj^-^^^l = u^^^kL- Then the ith entry of 
the vector v^-^ is less than or equal to the ith entry of the vector ujvi for alH G {1, . . . , m}, 
in other words VNi,i < UNi,i- We prove that 

'^Ni-kL+kp,i ^ '^Ni-kL+kp,ii 

by induction on p. The case p = provides the base case. Suppose the result holds for 
p < P. Then we have 

m rn 

VNi-kL+kP,i — ^N-i-kL+kP,ii'^^O.ijVNi-kL+k{P-l),j < 0'ijVNi-kL+k{P-l),j 

i=l j=l 

m 

— ^ 0'ijUNi-kL+k{P-l),j — UN-i-kL+kP,i, 

since bN^^kL+kP,ii G [0,1] and a^j > for all i.j G {l,...,m}. Thus we have shown 
in this case that VNi,i < UNi,i for all i G {1, . . . ,^1}. Notice that ^Jjv^^i = t/jVi,o,o,i and 
UNi,i = UNi,o,L,i- Further notice that this is true for any fixed L > and fixed Ni > kL. 
Thus Un,o,o,i < Unfi,L,i for alH G {1, . . . , m} and whenever L > and n > kL. So we 
have shown the case q — 0. 

Now we use the case g = as a base case in an induction argument on q. Suppose that 
Un,q,o,i ^ Un^q^L^i, then siucc all entries Ojj of our matrix A are non- negative and t^n,q,L is 
a vector with all non-negative entries, 

m m 
Un,q+l,0,i — y ] djj Un,q,0,j ^ 0'ijUn,q,L,j — t^n,q+l,L,i- 

i=i i=i 

For shorthand we define a function h : M"^ — )■ [0,oo) so that h{v) = {v,v). Fact 1 
gives us that there exists an G N so that h{A^v) < h{v) for all v G and all L > N. 
Thus fixing an L > iV and a g G N we have h{AiVn) < h^A^+Hn-kL) < h^AHn-kh)- 
Since each of the subsequences {h{A'^VnkL+a)}'^=i are monotone decreasing and bounded 
below by zero, they all converge. So lim„_!.oo h{A'^Vn) — h{A'^Vn-kL) — 0. By the squeeze 
theorem we get lim„^oo h{A'^Vn) — h^A^^'^Vn-kh) = 0. 

So the subsequences {h{vnkL+a)}^=i and {h{A^VnkL+a)}^=i are convergent and since 
lim„_j.oo h{vn) - h{A^Vn-kL) = we get 

lim h{VnkL+a) = ^a^ 1™ h{A^VnkL+a)- 

n— ^-oo n— >oo 

Now consider the sequence {vnkL+a}'^=i and let {vnjkL+a}j^i be a convergent subsequence 
with limj^.oo VrijkL+a = Wa- By what we have just shown it must be true that h{wa) — 
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h{A^Wa), but then by Fact 1 we have that Wa is in the span of the eigenvectors of A with 
corresponding eigenvalues whose absolute value is 1. Recall from the Perron- Frobenius 
theorem that there is only one such eigenvector and it is wi, the eigenvector associated 
to the eigenvalue 1. So Wa — cwi, where c is an arbitrary constant, and h{wa) = 2,a, also 

e [0, oo)"* as a consequence of our choice of initial conditions. Thus Wa — Wi ^ 



\Wl\ 



What this means is that the sequence {vnfcL+a}^i niust converge to Wa = Wi 
Suppose it does not, then for some e > there is a subsequence {vnakL+a}dLi so that 



'MdkL+a - wi [ jj-^ J 1 1 > e 



for all d E N. However we know that {vnakL+a}dLi is bounded and so it has a convergent 
subsequence. This means that {vnkL+a}'^!^=i has a convergent subsequence which does 

not converge to Wi (^jj^j- We have already shown that every convergent subsequence 
of {vnkL+a}'^=i couvcrges to Wi ■ Thus we have a contradiction. This proves that 

the sequence {vnkL+a}'^=i must converge to wi 



•u;i 

Now we want to apply Theorem 1 and Lemma 1. It is clear from the linear algebra 
discussion earlier that the sequence of matrices converges to a diagonal matrix 

with diagonal entries either 1 or 0. Since {-D"}^i converges, {A"}^^ converges. Thus 
j^n+i _ _^ Q p^Yso it is clear from this that {vn}'^=i is bounded and that there exists 
M > so that ll^^i^ll < Af||t;|| for all n G N and for all v E V. Now wc use Theorem 
1 to show that lim„^oo ||^n — ^'^n-fell = 0. So we apply Lemma 1 and we get that 

lim„^oo \ \Vn- Vn-k\ \ = 0. 

Combining the fact that {fnfcL+a}^! converges and lim„_^oo I l^n ~ Vn-k 1 1 = we get that 
if oi = 02 mod k then 

hm VnkL+ai = hm VnkL+a2- 
n—i-oo n—>oo 

Thus every solution must converge to a periodic solution of not necessarily prime period 
k. 

To construct a solution which is periodic with prime period k we use our eigenvector 
Wi associated with the eigenvalue 1 having strictly positive components. We choose 
initial conditions so that v„ = for n < 1 and n ^ 1 — k, and vi-k — wi. This is a 
periodic solution of prime period k. This concludes our proof. □ 

7. A Periodic Tetrachotomy Result 

Now we combine all of our work to give some preliminary examples of periodic tetra- 
chotomy behavior for systems of two rational difference equations. 

Theorem 4. Let k E N, k ^ 1 and consider the k*^ order system of two rational 
difference equations 

/3kXn-k + IkVn-k 



n e N, 
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Vn = fczri fcZi , n G 1^, 

with nonnegative parameters and nonnegative initial conditions. Define a matrix 

( Pk Ik 
\ 4 efc 

Assume 5k = Jk, then this system exhibits the following tetrachotomy behavior. 

I Suppose the spectral radius of A is less than 1, then every solution converges to the 
unique equilibrium. 

II Suppose the spectral radius of A is equal to 1 and —1 is not an eigenvalue of A, 
then every solution converges to a periodic solution of not necessarily prime period 
k. Furthermore in this case there exist periodic solutions with prime period k. 

III Suppose the spectral radius of A is equal to 1 and —1 is an eigenvalue of A, then 
every solution converges to a periodic solution of not necessarily prime period 2k. 
Furthermore in this case there exist periodic solutions with prime period 2k. 

IV Suppose the spectral radius of A is greater than 1 then there exist solutions where 
Xn + yn is unbounded. 

Proof. To begin we rewrite our system using matrix notation, as was done in Section 4. 
We let 

= r 1 , A 

Vn 



and 



where 









Ik 


4 









1 

















Our system then becomes 

Vn = BnAVn-k, U EN. 

Now case (I) follows directly from Theorem 1. Also case (IV) follows directly from 
Theorem 2. In case (II) the subcase where 7fc, 4, efc > is covered by Theorem 3. 
Thus all that remains to prove is case (II) in the subcase where one of the constants is 
and case (III). 

Recall that the solutions for the eigenvalues of a 2 x 2 matrix A can be written as 



A = i (^tr{A) ± Vtr2(A) - 4det{A) 



This computation is fairly straightforward; it appears as an exercise on page 39 in |15] . 
With our definition of A this becomes 



A = i + Cfc ± V(/3fe - efe)' + ^IkSk) . 
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we may write A in Jordan normal form and we get 

For the sake of notation we will sometimes use 

Ai 



D 



A 



Now consider a function : — )■ [0, oo) where h{v) = {v,v) the dot product of v with 
itself. Let wi and W2 be the unit eigenvectors of A, since A is diagonalizable wi and W2 
form an orthonormal basis for M^. It follows that for all e M^, 



2 2 

i=l i=l 
2 2 



^Y.\\,\%v,w,)\' <Y,\{^^^^)? = = h{v). 

1=1 i=l 

Also it is clear here that if h{Av) = h{v) and | Aj| < 1 for some i G {1, 2} then (f , Wi) = 0. 
In other words Fact 1 applies here, and since A is diagonalizable it applies for L > N — 1. 
Let Un = Avn-k, then we have 

l3kVn-k,l + lkVn-k,2 ^ n , jr 

"Vn.l = =;fcri < PkVn-k,l + lkVn-k,2 = ^71,1 

since a^, Vn-j G [0, oo)^. Also we have 

_ ^kVn-k.l + (kVn-k,2 ^ c- _ jj 

Vn,2 — S ('fc'Wn-fe.l + ^kVn-k,2 — l^n,2 

since qj,Vn-j G [0, oo)^. So < Un,i for all z G {1,2} where [/„ = Avn-k- Thus 
/i(-4v„) < < h{Avn-k) < h{vn-k)- 

Since each of the subsequences {h(ynk+a)}^=i 
and {h{Avnk+a)}'^=i are monotone decreasing and bounded below by 0, they all con- 
verge. So lim„_^oo h{vn) — h{vn-k) = 0. By the squeeze theorem we get lim„_^oo h{vn) — 
h{Avn-k) — 0. Thus we have 

lim h{vnk+a) ^ ^a= hm h{AVnk+a)- 

Now to finish case (II) our argument must break into several cases. 
Recall that in case (II) our eigenvalues must be real and not equal to —1 and must have 
the property that |Aj| < 1 since the spectral radius is 1 in this case. So we have two 
options, both Ai and A2 are equal to 1, or one of the two eigenvalues Aj has the property 
that |Aj| < 1. In the case where one of the two eigenvalues Aj has the property that 
I Ail < 1 we will use an argument almost identical to the argument which we applied in 
Theorem 3. If we are in this case one of the two eigenvalues, call it A2, has the property 
that IA2I < 1. If it is actually Ai with |Ai| < 1 then the argument is similar and will 
be omitted. Now consider the sequence {vnk+a}'^=i and let {vnjk+a}^=i be a convergent 
subsequence with limj^.oo ^n^k+a = ^a- By what we have just shown it must be true in 
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this case that h{wa) = h{Awa), but then by Fact 1 we have that Wa is in the span of the 
eigenvectors of A with corresponding eigenvalues whose absolute value is 1. Notice that 
in this case there is only one such eigenvector and it is wi, the eigenvector associated to 
the eigenvalue 1. So Wa — cwi, where c is an arbitrary constant, and h{wa) = S^a, also 

ifa e [0, oo)"^ as a consequence of our choice of initial conditions. Thus Wa — Wi ■ 

What this means is that the sequence {vnk+a}'^=i niust converge to Wa = Wi 
Suppose it does not, then for some e > there is a subsequence {vnak+a}dLi so that 



\Vndk+a - U)i { ) II > e 

for all c? e N. However we know that {fn^fc+ajdLi is bounded and so it has a convergent 
subsequence. This means that {vnk+a}'^=i has a convergent subsequence which does not 

converge to wi We have already shown that every convergent subsequence of 

{vnk+a}'^=i converges to Wi ^jy^j- Thus we have a contradiction. This proves that the 

sequence {vnk+a}'^=i must converge to Wi ^jj^^ ■ Thus in this case every solution must 
converge to a periodic solution of not necessarily prime period k. 

Now let us consider the other option in case (II). Assume that we have 5^ = 7^ = in 
case (II). Then we have forO<A<l, , or A — l)' focus 

on the recursive equations for x„ and y„, we get that 

Min-k ^ ^^T 

So we obtain the following recursive inequalities 

Xn < Xn-k, neN, 

Vn < ><yn-k, neN. 

So the subsequences {xnk+a}'^^=i ^■iid {ynk+a}'i^=i ^re all monotone decreasing and bounded 
below by zero, so they all converge and clearly — > 0. 
Or we have 

" — 1 ■sr^k—l r? v^fc— 1 ^ , c in, 

1 + l^j=l ^jXn-j + Z^j=l ^jVn-j 

Vn — J— i i — 1 , neN. 

So we obtain the following recursive inequalities 

Xn < Aa:„_fe, neN, 
Vn < Vn-k, neN. 
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So the subsequences {xnk+a}'^=i {ynk+a}'^=i are all monotone decreasing and bounded 
below by zero, so they all converge and clearly x„ — ?> 0. To construct a periodic solution 
take the initial conditions to be f„ = for n < 1 and n ^ 1 — k and vi^k = (0, 1) or 
vi-k = (1, 0) depending on the case. This is a periodic solution of prime period k. This 
concludes the proof of case (II). 

Now we will prove case (III). Now suppose /3k + ek > and —1 is an eigenvalue. Then 
we must have 



However since we have assumed that the spectral radius is 1 in this case that is impossi- 
ble. Thus /Sfc + efc < and we know from assumption that + > 0. Thus + = 
and in case (III) both —1 and 1 are eigenvalues. So in case (III) we have 

7fc 
^ 

So in case (III) we have the following system of rational difference equations 

IkVn-k ^ 

Vn = feZi feZi ; n e rsj. 

7fc(l + J2j=l DjXn-j + J2j=l EjVn-j) 

Thus we have the following recursive inequalities 

•^n — •^n—2ki 



Vn < yn-2k- 

So the subsequences {yn2k+a}'^=i and {a;„2fe+a}^i are all monotone decreasing and 
bounded below by zero, so they all converge. Thus we have shown that in case (III) 
every solution converges to a periodic solution of not necessarily prime period 2k. Since 
in case (III) we have 

choose initial conditions where = for n < 1 and n ^ 1 — k and 

vi-k = ( 6 ) ' 

where a,b & [0, oo) and a ^ 7^6 and the solution given by these initial conditions is a 
periodic solution of prime period 2k. This concludes the proof of case (III). 

□ 
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8. Conclusion 

We have created some analogues for trichotomy behavior for systems of rational differ- 
ence equations, but we have barely scratched the surface. There are literally thousands 
of special cases of systems of rational difference equations of order greater than one 
to explore. This paper leaves several questions for further study. Are there any other 
examples of periodic tetrachotomy behavior for systems of two rational difference equa- 
tions? Is it possible to make analogues to other trichotomy results in the literature? 
The main idea to take away from this article is that in some cases it is useful to reframe 
a problem about systems of rational difference equations as a problem about recursive 
systems on vector spaces. Doing this allows one to utilize the powerful tools available 
in linear algebra. 
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